matrix(96,40) estim=na
matrix(6,8) test_bs=na
matrix(4,8) test_asy=na

delete tb_*

'Hong Kong
'base
equation hk1.arch(0,1,m=2000,h) sqr(hk_hilo) @ hk_hilo(-1) 
equation hk1.arch(0,1,m=2000,h,b,s) sqr(hk_hilo)@ hk_hilo(-1) 
equation hk1.arch(0,1,m=2000,h,s) sqr(hk_hilo)@ hk_hilo(-1) 
estim(1,1)=hk1.@coefs(1)
estim(2,1)=hk1.@tstats(1)
estim(1,2)=hk1.@coefs(2)
estim(2,2)=hk1.@tstats(2)
estim(1,3)=hk1.@coefs(3)
estim(2,3)=hk1.@tstats(3)
estim(1,21)=hk1.@logl
estim(1,22)=hk1.@schwarz
estim(1,23)=hk1.@aic
hk1.makegarch mu_hk1
genr hk1_res=hk_hilo/mu_hk1
freeze(tb_1) hk1_res.ident(12)
estim(1,24)=tb_1(18,6)
estim(2,24)=tb_1(18,7)
freeze(tb_2) hk1_res^2.ident(12)
estim(1,25)=tb_2(18,6)
estim(2,25)=tb_2(18,7)
delete tb_*




' MOD EA (asymmetric effects during the crisis)

equation hk3a.arch(0,1,m=2000,h) sqr(hk_hilo) @ hk_hilo(-1) in_hilo(-1)  ko_hilo(-1) ma_hilo(-1)   ph_hilo(-1)  si_hilo(-1)   ta_hilo(-1) th_hilo(-1) hk_dpos(-1)*dc(-1) hk_dneg(-1)*dc(-1) 
'
equation hk3a.arch(0,1,m=2000,h,b,s) sqr(hk_hilo) @ hk_hilo(-1) in_hilo(-1)  ko_hilo(-1) ma_hilo(-1)   ph_hilo(-1)  si_hilo(-1)   ta_hilo(-1) th_hilo(-1) hk_dpos(-1)*dc(-1) hk_dneg(-1)*dc(-1) 
'
equation hk3a.arch(0,1,m=2000,h,s) sqr(hk_hilo) @ hk_hilo(-1) in_hilo(-1)  ko_hilo(-1) ma_hilo(-1)   ph_hilo(-1)  si_hilo(-1)   ta_hilo(-1) th_hilo(-1) hk_dpos(-1)*dc(-1) hk_dneg(-1)*dc(-1) 

estim(3,1)=hk3a.@coefs(1)
estim(4,1)=hk3a.@tstats(1)
estim(3,2)=hk3a.@coefs(2)
estim(4,2)=hk3a.@tstats(2)
estim(3,3)=hk3a.@coefs(3)
estim(4,3)=hk3a.@tstats(3)
estim(3,5)=hk3a.@coefs(4)
estim(4,5)=hk3a.@tstats(4)
estim(3,7)=hk3a.@coefs(5)
estim(4,7)=hk3a.@tstats(5)
estim(3,9)=hk3a.@coefs(6)
estim(4,9)=hk3a.@tstats(6)
estim(3,11)=hk3a.@coefs(7)
estim(4,11)=hk3a.@tstats(7)
estim(3,13)=hk3a.@coefs(8)
estim(4,13)=hk3a.@tstats(8)
estim(3,15)=hk3a.@coefs(9)
estim(4,15)=hk3a.@tstats(9)
estim(3,17)=hk3a.@coefs(10)
estim(4,17)=hk3a.@tstats(10)
estim(3,29)=hk3a.@coefs(11)
estim(4,29)=hk3a.@tstats(11)
estim(3,30)=hk3a.@coefs(12)
estim(4,30)=hk3a.@tstats(12)
estim(3,21)=hk3a.@logl
estim(3,22)=hk3a.@schwarz
estim(3,23)=hk3a.@aic
hk3a.makegarch mu_hk3a
genr hk3a_res=hk_hilo/mu_hk3a
freeze(tb_1) hk3a_res.ident(12)
estim(3,24)=tb_1(18,6)
estim(4,24)=tb_1(18,7)
freeze(tb_2) hk3a_res^2.ident(12)
estim(3,25)=tb_2(18,6)
estim(4,25)=tb_2(18,7)
delete tb_*



'Indonesia
'base
equation in1.arch(0,1,m=2000,h) sqr(in_hilo) @ in_hilo(-1)
equation in1.arch(0,1,m=2000,h,b,s) sqr(in_hilo)@ in_hilo(-1)
equation in1.arch(0,1,m=2000,h,s) sqr(in_hilo)@ in_hilo(-1)
estim(5,1)=in1.@coefs(1)
estim(6,1)=in1.@tstats(1)
estim(5,2)=in1.@coefs(2)
estim(6,2)=in1.@tstats(2)
estim(5,5)=in1.@coefs(3)
estim(6,5)=in1.@tstats(3)
estim(5,21)=in1.@logl
estim(5,22)=in1.@schwarz
estim(5,23)=in1.@aic
in1.makegarch mu_in1
genr in1_res=in_hilo/mu_in1
freeze(tb_1) in1_res.ident(12)
estim(5,24)=tb_1(18,6)
estim(6,24)=tb_1(18,7)
freeze(tb_2) in1_res^2.ident(12)
estim(5,25)=tb_2(18,6)
estim(6,25)=tb_2(18,7)
delete tb_*

'selected: completo con dc, dc1 e  interazione mkt*dc (MOD EXD)

equation in4e.arch(0,1,m=2000,h) sqr(in_hilo) @ hk_hilo(-1) hk_hilo(-1)*dc(-1) in_hilo(-1) in_hilo(-1)*dc(-1) ko_hilo(-1) ko_hilo(-1)*dc(-1) ma_hilo(-1) ma_hilo(-1)*dc(-1) ph_hilo(-1) ph_hilo(-1)*dc(-1)  si_hilo(-1) si_hilo(-1)*dc(-1) ta_hilo(-1) ta_hilo(-1)*dc(-1)  th_hilo(-1) th_hilo(-1)*dc(-1) dc(-1) dc1(-1) 

equation in4e.arch(0,1,m=2000,h,b,s) sqr(in_hilo)  @ hk_hilo(-1) hk_hilo(-1)*dc(-1) in_hilo(-1) in_hilo(-1)*dc(-1) ko_hilo(-1) ko_hilo(-1)*dc(-1) ma_hilo(-1) ma_hilo(-1)*dc(-1) ph_hilo(-1) ph_hilo(-1)*dc(-1)  si_hilo(-1) si_hilo(-1)*dc(-1) ta_hilo(-1) ta_hilo(-1)*dc(-1)  th_hilo(-1) th_hilo(-1)*dc(-1) dc(-1) dc1(-1) 


equation in4e.arch(0,1,m=2000,h,s) sqr(in_hilo)  @ hk_hilo(-1) hk_hilo(-1)*dc(-1) in_hilo(-1) in_hilo(-1)*dc(-1) ko_hilo(-1) ko_hilo(-1)*dc(-1) ma_hilo(-1) ma_hilo(-1)*dc(-1) ph_hilo(-1) ph_hilo(-1)*dc(-1)  si_hilo(-1) si_hilo(-1)*dc(-1) ta_hilo(-1) ta_hilo(-1)*dc(-1)  th_hilo(-1) th_hilo(-1)*dc(-1) dc(-1) dc1(-1) 
estim(7,1)=in4e.@coefs(1)
estim(8,1)=in4e.@tstats(1)
estim(7,2)=in4e.@coefs(2)
estim(8,2)=in4e.@tstats(2)
estim(7,3)=in4e.@coefs(3)
estim(8,3)=in4e.@tstats(3)
estim(7,4)=in4e.@coefs(4)
estim(8,4)=in4e.@tstats(4)
estim(7,5)=in4e.@coefs(5)
estim(8,5)=in4e.@tstats(5)
estim(7,6)=in4e.@coefs(6)
estim(8,6)=in4e.@tstats(6)
estim(7,7)=in4e.@coefs(7)
estim(8,7)=in4e.@tstats(7)
estim(7,8)=in4e.@coefs(8)
estim(8,8)=in4e.@tstats(8)
estim(7,9)=in4e.@coefs(9)
estim(8,9)=in4e.@tstats(9)
estim(7,10)=in4e.@coefs(10)
estim(8,10)=in4e.@tstats(10)
estim(7,11)=in4e.@coefs(11)
estim(8,11)=in4e.@tstats(11)
estim(7,12)=in4e.@coefs(12)
estim(8,12)=in4e.@tstats(12)
estim(7,13)=in4e.@coefs(13)
estim(8,13)=in4e.@tstats(13)
estim(7,14)=in4e.@coefs(14)
estim(8,14)=in4e.@tstats(14)
estim(7,15)=in4e.@coefs(15)
estim(8,15)=in4e.@tstats(15)
estim(7,16)=in4e.@coefs(16)
estim(8,16)=in4e.@tstats(16)
estim(7,17)=in4e.@coefs(17)
estim(8,17)=in4e.@tstats(17)
estim(7,18)=in4e.@coefs(18)
estim(8,18)=in4e.@tstats(18)
estim(7,19)=in4e.@coefs(19)
estim(8,19)=in4e.@tstats(19)
estim(7,20)=in4e.@coefs(20)
estim(8,20)=in4e.@tstats(20)
estim(7,21)=in4e.@logl
estim(7,22)=in4e.@schwarz
estim(7,23)=in4e.@aic
in4e.makegarch mu_in4e
genr in4e_res=in_hilo/mu_in4e
freeze(tb_1) in4e_res.ident(12)
estim(7,24)=tb_1(18,6)
estim(8,24)=tb_1(18,7)
freeze(tb_2) in4e_res^2.ident(12)
estim(7,25)=tb_2(18,6)
estim(8,25)=tb_2(18,7)
delete tb_*





'Korea
'base
equation ko1.arch(0,1,m=2000,h) sqr(ko_hilo) @ ko_dpos(-1) ko_dneg(-1)
equation ko1.arch(0,1,m=2000,h,b,s) sqr(ko_hilo)@ ko_dpos(-1) ko_dneg(-1)
equation ko1.arch(0,1,m=2000,h,s) sqr(ko_hilo)@ ko_dpos(-1) ko_dneg(-1)
estim(9,1)=ko1.@coefs(1)
estim(10,1)=ko1.@tstats(1)
estim(9,2)=ko1.@coefs(2)
estim(10,2)=ko1.@tstats(2)
estim(9,7)=ko1.@coefs(3)
estim(10,7)=ko1.@tstats(3)
estim(9,31)=ko1.@coefs(4)
estim(10,31)=ko1.@tstats(4)
estim(9,21)=ko1.@logl
estim(9,22)=ko1.@schwarz
estim(9,23)=ko1.@aic
ko1.makegarch mu_ko1
genr ko1_res=ko_hilo/mu_ko1
freeze(tb_1) ko1_res.ident(12)
estim(9,24)=tb_1(18,6)
estim(10,24)=tb_1(18,7)
freeze(tb_2) ko1_res^2.ident(12)
estim(9,25)=tb_2(18,6)
estim(10,25)=tb_2(18,7)
delete tb_*




'completo con dc, dc1,  interazione mkt*dc, effetti asimmetrici fuori dalla crisi (MOD EXDA)

equation ko4ea.arch(0,1,m=2000,h) sqr(ko_hilo) @ hk_hilo(-1) hk_hilo(-1)*dc(-1) in_hilo(-1) in_hilo(-1)*dc(-1) ko_dpos(-1) ko_hilo(-1)*dc(-1) ma_hilo(-1) ma_hilo(-1)*dc(-1) ph_hilo(-1) ph_hilo(-1)*dc(-1)  si_hilo(-1) si_hilo(-1)*dc(-1) ta_hilo(-1) ta_hilo(-1)*dc(-1)  th_hilo(-1) th_hilo(-1)*dc(-1) dc(-1) dc1(-1) ko_dneg(-1)  


equation ko4ea.arch(0,1,m=2000,h,b,s) sqr(ko_hilo)  @ hk_hilo(-1) hk_hilo(-1)*dc(-1) in_hilo(-1) in_hilo(-1)*dc(-1) ko_dpos(-1) ko_hilo(-1)*dc(-1) ma_hilo(-1) ma_hilo(-1)*dc(-1) ph_hilo(-1) ph_hilo(-1)*dc(-1)  si_hilo(-1) si_hilo(-1)*dc(-1) ta_hilo(-1) ta_hilo(-1)*dc(-1)  th_hilo(-1) th_hilo(-1)*dc(-1) dc(-1) dc1(-1) ko_dneg(-1)  


equation ko4ea.arch(0,1,m=2000,h,s) sqr(ko_hilo)  @ hk_hilo(-1) hk_hilo(-1)*dc(-1) in_hilo(-1) in_hilo(-1)*dc(-1) ko_dpos(-1) ko_hilo(-1)*dc(-1) ma_hilo(-1) ma_hilo(-1)*dc(-1) ph_hilo(-1) ph_hilo(-1)*dc(-1)  si_hilo(-1) si_hilo(-1)*dc(-1) ta_hilo(-1) ta_hilo(-1)*dc(-1)  th_hilo(-1) th_hilo(-1)*dc(-1) dc(-1) dc1(-1) ko_dneg(-1)  

estim(11,1)=ko4ea.@coefs(1)
estim(12,1)=ko4ea.@tstats(1)
estim(11,2)=ko4ea.@coefs(2)
estim(12,2)=ko4ea.@tstats(2)
estim(11,3)=ko4ea.@coefs(3)
estim(12,3)=ko4ea.@tstats(3)
estim(11,4)=ko4ea.@coefs(4)
estim(12,4)=ko4ea.@tstats(4)
estim(11,5)=ko4ea.@coefs(5)
estim(12,5)=ko4ea.@tstats(5)
estim(11,6)=ko4ea.@coefs(6)
estim(12,6)=ko4ea.@tstats(6)
estim(11,7)=ko4ea.@coefs(7)
estim(12,7)=ko4ea.@tstats(7)
estim(11,8)=ko4ea.@coefs(8)
estim(12,8)=ko4ea.@tstats(8)
estim(11,9)=ko4ea.@coefs(9)
estim(12,9)=ko4ea.@tstats(9)
estim(11,10)=ko4ea.@coefs(10)
estim(12,10)=ko4ea.@tstats(10)
estim(11,11)=ko4ea.@coefs(11)
estim(12,11)=ko4ea.@tstats(11)
estim(11,12)=ko4ea.@coefs(12)
estim(12,12)=ko4ea.@tstats(12)
estim(11,13)=ko4ea.@coefs(13)
estim(12,13)=ko4ea.@tstats(13)
estim(11,14)=ko4ea.@coefs(14)
estim(12,14)=ko4ea.@tstats(14)
estim(11,15)=ko4ea.@coefs(15)
estim(12,15)=ko4ea.@tstats(15)
estim(11,16)=ko4ea.@coefs(16)
estim(12,16)=ko4ea.@tstats(16)
estim(11,17)=ko4ea.@coefs(17)
estim(12,17)=ko4ea.@tstats(17)
estim(11,18)=ko4ea.@coefs(18)
estim(12,18)=ko4ea.@tstats(18)
estim(11,19)=ko4ea.@coefs(19)
estim(12,19)=ko4ea.@tstats(19)
estim(11,20)=ko4ea.@coefs(20)
estim(12,20)=ko4ea.@tstats(20)
estim(11,31)=ko4ea.@coefs(21)
estim(12,31)=ko4ea.@tstats(21)
estim(11,21)=ko4ea.@logl
estim(11,22)=ko4ea.@schwarz
estim(11,23)=ko4ea.@aic
ko4ea.makegarch mu_ko4ea
genr ko4ea_res=ko_hilo/mu_ko4ea
freeze(tb_1) ko4ea_res.ident(12)
estim(11,24)=tb_1(18,6)
estim(12,24)=tb_1(18,7)
freeze(tb_2) ko4ea_res^2.ident(12)
estim(11,25)=tb_2(18,6)
estim(12,25)=tb_2(18,7)
delete tb_*




'Malaysia

'base
equation ma1.arch(0,1,m=2000,h) sqr(ma_hilo) @ ma_hilo(-1) ma_hilo(-2) 
equation ma1.arch(0,1,m=2000,h,b,s) sqr(ma_hilo)@ ma_hilo(-1) ma_hilo(-2) 
equation ma1.arch(0,1,m=2000,h,s) sqr(ma_hilo)@ ma_hilo(-1) ma_hilo(-2) 
estim(13,1)=ma1.@coefs(1)
estim(14,1)=ma1.@tstats(1)
estim(13,2)=ma1.@coefs(2)
estim(14,2)=ma1.@tstats(2)
estim(13,9)=ma1.@coefs(3)
estim(14,9)=ma1.@tstats(3)
estim(13,26)=ma1.@coefs(4)
estim(14,26)=ma1.@tstats(4)
estim(13,21)=ma1.@logl
estim(13,22)=ma1.@schwarz
estim(13,23)=ma1.@aic
ma1.makegarch mu_ma1
genr ma1_res=ma_hilo/mu_ma1
freeze(tb_1) ma1_res.ident(12)
estim(13,24)=tb_1(18,6)
estim(14,24)=tb_1(18,7)
freeze(tb_2) ma1_res^2.ident(12)
estim(13,25)=tb_2(18,6)
estim(14,25)=tb_2(18,7)
delete tb_*




' completo con dc(-1) senza interazione  (MOD EDC)
equation ma5.arch(0,1,m=2000,h) sqr(ma_hilo) @ hk_hilo(-1) in_hilo(-1)    ko_hilo(-1) ma_hilo(-1) ph_hilo(-1)  si_hilo(-1)  ta_hilo(-1)  th_hilo(-1) dc(-1)  ma_hilo(-2) 
'
equation ma5.arch(0,1,m=2000,h,b,s) sqr(ma_hilo) @  hk_hilo(-1) in_hilo(-1)    ko_hilo(-1) ma_hilo(-1) ph_hilo(-1)  si_hilo(-1)  ta_hilo(-1)  th_hilo(-1) dc(-1)  ma_hilo(-2) 
'
equation ma5.arch(0,1,m=2000,h,s) sqr(ma_hilo) @ hk_hilo(-1) in_hilo(-1)    ko_hilo(-1) ma_hilo(-1) ph_hilo(-1)  si_hilo(-1)  ta_hilo(-1)  th_hilo(-1) dc(-1)  ma_hilo(-2) 
estim(15,1)=ma5.@coefs(1)
estim(16,1)=ma5.@tstats(1)
estim(15,2)=ma5.@coefs(2)
estim(16,2)=ma5.@tstats(2)
estim(15,3)=ma5.@coefs(3)
estim(16,3)=ma5.@tstats(3)
estim(15,5)=ma5.@coefs(4)
estim(16,5)=ma5.@tstats(4)
estim(15,7)=ma5.@coefs(5)
estim(16,7)=ma5.@tstats(5)
estim(15,9)=ma5.@coefs(6)
estim(16,9)=ma5.@tstats(6)
estim(15,11)=ma5.@coefs(7)
estim(16,11)=ma5.@tstats(7)
estim(15,13)=ma5.@coefs(8)
estim(16,13)=ma5.@tstats(8)
estim(15,15)=ma5.@coefs(9)
estim(16,15)=ma5.@tstats(9)
estim(15,17)=ma5.@coefs(10)
estim(16,17)=ma5.@tstats(10)
estim(15,19)=ma5.@coefs(11)
estim(16,19)=ma5.@tstats(11)
estim(15,26)=ma5.@coefs(12)
estim(16,26)=ma5.@tstats(12)
estim(15,21)=ma5.@logl
estim(15,22)=ma5.@schwarz
estim(15,23)=ma5.@aic
ma5.makegarch mu_ma5
genr ma5_res=ma_hilo/mu_ma5
freeze(tb_1) ma5_res.ident(12)
estim(15,24)=tb_1(18,6)
estim(16,24)=tb_1(18,7)
freeze(tb_2) ma5_res^2.ident(12)
estim(15,25)=tb_2(18,6)
estim(16,25)=tb_2(18,7)
delete tb_*





'Philippines
'base
equation ph1.arch(0,1,m=2000,h) sqr(ph_hilo) @ ph_hilo(-1) 
equation ph1.arch(0,1,m=2000,h,b,s) sqr(ph_hilo)@ ph_hilo(-1) 
equation ph1.arch(0,1,m=2000,h,s) sqr(ph_hilo)@ ph_hilo(-1) 
estim(17,1)=ph1.@coefs(1)
estim(18,1)=ph1.@tstats(1)
estim(17,2)=ph1.@coefs(2)
estim(18,2)=ph1.@tstats(2)
estim(17,11)=ph1.@coefs(3)
estim(18,11)=ph1.@tstats(3)
estim(17,21)=ph1.@logl
estim(17,22)=ph1.@schwarz
estim(17,23)=ph1.@aic
ph1.makegarch mu_ph1
genr ph1_res=ph_hilo/mu_ph1
freeze(tb_1) ph1_res.ident(12)
estim(17,24)=tb_1(18,6)
estim(18,24)=tb_1(18,7)
freeze(tb_2) ph1_res^2.ident(12)
estim(17,25)=tb_2(18,6)
estim(18,25)=tb_2(18,7)
delete tb_*



' completo con dc(-1) senza interazione (MOD EDC)

equation ph5.arch(0,1,m=2000,h) sqr(ph_hilo) @  hk_hilo(-1) in_hilo(-1)   ko_hilo(-1) ma_hilo(-1)  ph_hilo(-1)  si_hilo(-1)   ta_hilo(-1) th_hilo(-1) dc(-1) 
'
equation ph5.arch(0,1,m=2000,h,b,s) sqr(ph_hilo) @ hk_hilo(-1) in_hilo(-1)   ko_hilo(-1) ma_hilo(-1)  ph_hilo(-1)  si_hilo(-1)   ta_hilo(-1) th_hilo(-1) dc(-1) 
'
equation ph5.arch(0,1,m=2000,h,s) sqr(ph_hilo) @ hk_hilo(-1) in_hilo(-1)   ko_hilo(-1) ma_hilo(-1)  ph_hilo(-1)  si_hilo(-1)   ta_hilo(-1) th_hilo(-1) dc(-1) 
estim(19,1)=ph5.@coefs(1)
estim(20,1)=ph5.@tstats(1)
estim(19,2)=ph5.@coefs(2)
estim(20,2)=ph5.@tstats(2)
estim(19,3)=ph5.@coefs(3)
estim(20,3)=ph5.@tstats(3)
estim(19,5)=ph5.@coefs(4)
estim(20,5)=ph5.@tstats(4)
estim(19,7)=ph5.@coefs(5)
estim(20,7)=ph5.@tstats(5)
estim(19,9)=ph5.@coefs(6)
estim(20,9)=ph5.@tstats(6)
estim(19,11)=ph5.@coefs(7)
estim(20,11)=ph5.@tstats(7)
estim(19,13)=ph5.@coefs(8)
estim(20,13)=ph5.@tstats(8)
estim(19,15)=ph5.@coefs(9)
estim(20,15)=ph5.@tstats(9)
estim(19,17)=ph5.@coefs(10)
estim(20,17)=ph5.@tstats(10)
estim(19,19)=ph5.@coefs(11)
estim(20,19)=ph5.@tstats(11)
estim(19,21)=ph5.@logl
estim(19,22)=ph5.@schwarz
estim(19,23)=ph5.@aic
ph5.makegarch mu_ph5
genr ph5_res=ph_hilo/mu_ph5
freeze(tb_1) ph5_res.ident(12)
estim(19,24)=tb_1(18,6)
estim(20,24)=tb_1(18,7)
freeze(tb_2) ph5_res^2.ident(12)
estim(19,25)=tb_2(18,6)
estim(20,25)=tb_2(18,7)
delete tb_*




'Singapore
'base
equation si1.arch(0,1,m=2000,h) sqr(si_hilo) @ si_hilo(-1) si_hilo(-2) 
equation si1.arch(0,1,m=2000,h,b,s) sqr(si_hilo)@ si_hilo(-1) si_hilo(-2) 
equation si1.arch(0,1,m=2000,h,s) sqr(si_hilo)@ si_hilo(-1) si_hilo(-2) 
estim(21,1)=si1.@coefs(1)
estim(22,1)=si1.@tstats(1)
estim(21,2)=si1.@coefs(2)
estim(22,2)=si1.@tstats(2)
estim(21,13)=si1.@coefs(3)
estim(22,13)=si1.@tstats(3)
estim(21,27)=si1.@coefs(4)
estim(22,27)=si1.@tstats(4)
estim(21,21)=si1.@logl
estim(21,22)=si1.@schwarz
estim(21,23)=si1.@aic
si1.makegarch mu_si1
genr si1_res=si_hilo/mu_si1
freeze(tb_1) si1_res.ident(12)
estim(21,24)=tb_1(18,6)
estim(22,24)=tb_1(18,7)
freeze(tb_2) si1_res^2.ident(12)
estim(21,25)=tb_2(18,6)
estim(22,25)=tb_2(18,7)
delete tb_*




' completo con effetti asimmetrici durante la crisi (MOD EA)
equation si3a.arch(0,1,m=2000,h) sqr(si_hilo) @ hk_hilo(-1) in_hilo(-1)  ko_hilo(-1) ma_hilo(-1) ph_hilo(-1)  si_hilo(-1)  ta_hilo(-1) th_hilo(-1) si_hilo(-2) si_dpos(-1)*dc(-1) si_dneg(-1)*dc(-1)
'
equation si3a.arch(0,1,m=2000,h,b,s) sqr(si_hilo) @ hk_hilo(-1) in_hilo(-1)  ko_hilo(-1) ma_hilo(-1) ph_hilo(-1)  si_hilo(-1)  ta_hilo(-1) th_hilo(-1) si_hilo(-2) si_dpos(-1)*dc(-1) si_dneg(-1)*dc(-1)
'
equation si3a.arch(0,1,m=2000,h,s) sqr(si_hilo) @ hk_hilo(-1) in_hilo(-1)  ko_hilo(-1) ma_hilo(-1) ph_hilo(-1)  si_hilo(-1)  ta_hilo(-1) th_hilo(-1) si_hilo(-2) si_dpos(-1)*dc(-1) si_dneg(-1)*dc(-1)

estim(23,1)=si3a.@coefs(1)
estim(24,1)=si3a.@tstats(1)
estim(23,2)=si3a.@coefs(2)
estim(24,2)=si3a.@tstats(2)
estim(23,3)=si3a.@coefs(3)
estim(24,3)=si3a.@tstats(3)
estim(23,5)=si3a.@coefs(4)
estim(24,5)=si3a.@tstats(4)
estim(23,7)=si3a.@coefs(5)
estim(24,7)=si3a.@tstats(5)
estim(23,9)=si3a.@coefs(6)
estim(24,9)=si3a.@tstats(6)
estim(23,11)=si3a.@coefs(7)
estim(24,11)=si3a.@tstats(7)
estim(23,13)=si3a.@coefs(8)
estim(24,13)=si3a.@tstats(8)
estim(23,15)=si3a.@coefs(9)
estim(24,15)=si3a.@tstats(9)
estim(23,17)=si3a.@coefs(10)
estim(24,17)=si3a.@tstats(10)
estim(23,27)=si3a.@coefs(11)
estim(24,27)=si3a.@tstats(11)
estim(23,32)=si3a.@coefs(12)
estim(24,32)=si3a.@tstats(12)
estim(23,33)=si3a.@coefs(13)
estim(24,33)=si3a.@tstats(13)
estim(23,21)=si3a.@logl
estim(23,22)=si3a.@schwarz
estim(23,23)=si3a.@aic
si3a.makegarch mu_si3a
genr si3a_res=si_hilo/mu_si3a
freeze(tb_1) si3a_res.ident(12)
estim(23,24)=tb_1(18,6)
estim(24,24)=tb_1(18,7)
freeze(tb_2) si3a_res^2.ident(12)
estim(23,25)=tb_2(18,6)
estim(24,25)=tb_2(18,7)
delete tb_*







'Taiwan
'base
equation ta1.arch(0,1,m=2000,h) sqr(ta_hilo) @ ta_dpos(-1) ta_dneg(-1)
equation ta1.arch(0,1,m=2000,h,b,s) sqr(ta_hilo)@ ta_dpos(-1) ta_dneg(-1)
equation ta1.arch(0,1,m=2000,h,s) sqr(ta_hilo)@ ta_dpos(-1) ta_dneg(-1)
estim(25,1)=ta1.@coefs(1)
estim(26,1)=ta1.@tstats(1)
estim(25,2)=ta1.@coefs(2)
estim(26,2)=ta1.@tstats(2)
estim(25,15)=ta1.@coefs(3)
estim(26,15)=ta1.@tstats(3)
estim(25,34)=ta1.@coefs(4)
estim(26,34)=ta1.@tstats(4)
estim(25,21)=ta1.@logl
estim(25,22)=ta1.@schwarz
estim(25,23)=ta1.@aic
ta1.makegarch mu_ta1
genr ta1_res=ta_hilo/mu_ta1
freeze(tb_1) ta1_res.ident(12)
estim(25,24)=tb_1(18,6)
estim(26,24)=tb_1(18,7)
freeze(tb_2) ta1_res^2.ident(12)
estim(25,25)=tb_2(18,6)
estim(26,25)=tb_2(18,7)
delete tb_*



'completo con interazione mkt*dc(-1) ed effetti asimmetrici su tutto il periodo e durante la crisi (MOD EXAA)

equation ta4aa.arch(0,1,m=2000,h) sqr(ta_hilo) @ hk_hilo(-1) hk_hilo(-1)*dc(-1) in_hilo(-1) in_hilo(-1)*dc(-1) ko_hilo(-1) ko_hilo(-1)*dc(-1) ma_hilo(-1) ma_hilo(-1)*dc(-1)  ph_hilo(-1) ph_hilo(-1)*dc(-1)  si_hilo(-1) si_hilo(-1)*dc(-1)  ta_dpos(-1) ta_dpos(-1)*dc(-1) th_hilo(-1) th_hilo(-1)*dc(-1) ta_dneg(-1) ta_dneg(-1)*dc(-1)
'
equation ta4aa.arch(0,1,m=2000,h,b,s) sqr(ta_hilo) @ hk_hilo(-1) hk_hilo(-1)*dc(-1) in_hilo(-1) in_hilo(-1)*dc(-1)  ko_hilo(-1) ko_hilo(-1)*dc(-1) ma_hilo(-1) ma_hilo(-1)*dc(-1) ph_hilo(-1) ph_hilo(-1)*dc(-1)  si_hilo(-1) si_hilo(-1)*dc(-1) ta_dpos(-1) ta_dpos(-1)*dc(-1) th_hilo(-1) th_hilo(-1)*dc(-1) ta_dneg(-1) ta_dneg(-1)*dc(-1)
'
equation ta4aa.arch(0,1,m=2000,h,s) sqr(ta_hilo) @ hk_hilo(-1) hk_hilo(-1)*dc(-1) in_hilo(-1) in_hilo(-1)*dc(-1) ko_hilo(-1) ko_hilo(-1)*dc(-1) ma_hilo(-1) ma_hilo(-1)*dc(-1) ph_hilo(-1) ph_hilo(-1)*dc(-1)  si_hilo(-1) si_hilo(-1)*dc(-1)  ta_dpos(-1) ta_dpos(-1)*dc(-1) th_hilo(-1) th_hilo(-1)*dc(-1) ta_dneg(-1) ta_dneg(-1)*dc(-1)
estim(27,1)=ta4aa.@coefs(1)
estim(28,1)=ta4aa.@tstats(1)
estim(27,2)=ta4aa.@coefs(2)
estim(28,2)=ta4aa.@tstats(2)
estim(27,3)=ta4aa.@coefs(3)
estim(28,3)=ta4aa.@tstats(3)
estim(27,4)=ta4aa.@coefs(4)
estim(28,4)=ta4aa.@tstats(4)
estim(27,5)=ta4aa.@coefs(5)
estim(28,5)=ta4aa.@tstats(5)
estim(27,6)=ta4aa.@coefs(6)
estim(28,6)=ta4aa.@tstats(6)
estim(27,7)=ta4aa.@coefs(7)
estim(28,7)=ta4aa.@tstats(7)
estim(27,8)=ta4aa.@coefs(8)
estim(28,8)=ta4aa.@tstats(8)
estim(27,9)=ta4aa.@coefs(9)
estim(28,9)=ta4aa.@tstats(9)
estim(27,10)=ta4aa.@coefs(10)
estim(28,10)=ta4aa.@tstats(10)
estim(27,11)=ta4aa.@coefs(11)
estim(28,11)=ta4aa.@tstats(11)
estim(27,12)=ta4aa.@coefs(12)
estim(28,12)=ta4aa.@tstats(12)
estim(27,13)=ta4aa.@coefs(13)
estim(28,13)=ta4aa.@tstats(13)
estim(27,14)=ta4aa.@coefs(14)
estim(28,14)=ta4aa.@tstats(14)
estim(27,15)=ta4aa.@coefs(15)
estim(28,15)=ta4aa.@tstats(15)
estim(27,16)=ta4aa.@coefs(16)
estim(28,16)=ta4aa.@tstats(16)
estim(27,17)=ta4aa.@coefs(17)
estim(28,17)=ta4aa.@tstats(17)
estim(27,18)=ta4aa.@coefs(18)
estim(28,18)=ta4aa.@tstats(18)
estim(27,34)=ta4aa.@coefs(19)
estim(28,34)=ta4aa.@tstats(19)
estim(27,35)=ta4aa.@coefs(20)
estim(28,35)=ta4aa.@tstats(20)
estim(27,21)=ta4aa.@logl
estim(27,22)=ta4aa.@schwarz
estim(27,23)=ta4aa.@aic
ta4aa.makegarch mu_ta4aa
genr ta4aa_res=ta_hilo/mu_ta4aa
freeze(tb_1) ta4aa_res.ident(12)
estim(27,24)=tb_1(18,6)
estim(28,24)=tb_1(18,7)
freeze(tb_2) ta4aa_res^2.ident(12)
estim(27,25)=tb_2(18,6)
estim(28,25)=tb_2(18,7)
delete tb_*


'Thailand
'base
equation th1.arch(0,1,m=2000,h) sqr(th_hilo) @ th_hilo(-1) th_hilo(-2) 
equation th1.arch(0,1,m=2000,h,b,s) sqr(th_hilo)@ th_hilo(-1) th_hilo(-2) 
equation th1.arch(0,1,m=2000,h,s) sqr(th_hilo)@ th_hilo(-1) th_hilo(-2) 
estim(29,1)=th1.@coefs(1)
estim(30,1)=th1.@tstats(1)
estim(29,2)=th1.@coefs(2)
estim(30,2)=th1.@tstats(2)
estim(29,17)=th1.@coefs(3)
estim(30,17)=th1.@tstats(3)
estim(29,28)=th1.@coefs(4)
estim(30,28)=th1.@tstats(4)
estim(29,21)=th1.@logl
estim(29,22)=th1.@schwarz
estim(29,23)=th1.@aic
th1.makegarch mu_th1
genr th1_res=th_hilo/mu_th1
freeze(tb_1) th1_res.ident(12)
estim(29,24)=tb_1(18,6)
estim(30,24)=tb_1(18,7)
freeze(tb_2) th1_res^2.ident(12)
estim(29,25)=tb_2(18,6)
estim(30,25)=tb_2(18,7)
delete tb_*



' completo con effetti asimmetrici durante la crisi (MOD EA)
equation th3a.arch(0,1,m=2000,h) sqr(th_hilo) @ hk_hilo(-1) in_hilo(-1)  ko_hilo(-1)  ma_hilo(-1)  ph_hilo(-1)  si_hilo(-1)  ta_hilo(-1) th_hilo(-1) th_hilo(-2) th_dpos(-1)*dc(-1) th_dneg(-1)*dc(-1)
'
equation th3a.arch(0,1,m=2000,h,b,s) sqr(th_hilo) @ hk_hilo(-1) in_hilo(-1)  ko_hilo(-1)  ma_hilo(-1)  ph_hilo(-1)  si_hilo(-1)  ta_hilo(-1) th_hilo(-1) th_hilo(-2) th_dpos(-1)*dc(-1) th_dneg(-1)*dc(-1)
'
equation th3a.arch(0,1,m=2000,h,s) sqr(th_hilo) @ hk_hilo(-1) in_hilo(-1)  ko_hilo(-1)  ma_hilo(-1)  ph_hilo(-1)  si_hilo(-1)  ta_hilo(-1) th_hilo(-1) th_hilo(-2) th_dpos(-1)*dc(-1) th_dneg(-1)*dc(-1)

estim(31,1)=th3a.@coefs(1)
estim(32,1)=th3a.@tstats(1)
estim(31,2)=th3a.@coefs(2)
estim(32,2)=th3a.@tstats(2)
estim(31,3)=th3a.@coefs(3)
estim(32,3)=th3a.@tstats(3)
estim(31,5)=th3a.@coefs(4)
estim(32,5)=th3a.@tstats(4)
estim(31,7)=th3a.@coefs(5)
estim(32,7)=th3a.@tstats(5)
estim(31,9)=th3a.@coefs(6)
estim(32,9)=th3a.@tstats(6)
estim(31,11)=th3a.@coefs(7)
estim(32,11)=th3a.@tstats(7)
estim(31,13)=th3a.@coefs(8)
estim(32,13)=th3a.@tstats(8)
estim(31,15)=th3a.@coefs(9)
estim(32,15)=th3a.@tstats(9)
estim(31,17)=th3a.@coefs(10)
estim(32,17)=th3a.@tstats(10)
estim(31,28)=th3a.@coefs(11)
estim(32,28)=th3a.@tstats(11)
estim(31,36)=th3a.@coefs(12)
estim(32,36)=th3a.@tstats(12)
estim(31,37)=th3a.@coefs(13)
estim(32,37)=th3a.@tstats(13)
estim(31,21)=th3a.@logl
estim(31,22)=th3a.@schwarz
estim(31,23)=th3a.@aic
th3a.makegarch mu_th3a
genr th3a_res=th_hilo/mu_th3a
freeze(tb_1) th3a_res.ident(12)
estim(31,24)=tb_1(18,6)
estim(32,24)=tb_1(18,7)
freeze(tb_2) th3a_res^2.ident(12)
estim(31,25)=tb_2(18,6)
estim(32,25)=tb_2(18,7)


group eps hk3a_res in4e_res ko4ea_res ma5_res ph5_res si3a_res ta4aa_res th3a_res

group mu mu_hk3a mu_in4e mu_ko4ea mu_ma5 mu_ph5 mu_si3a mu_ta4aa mu_th3a


delete tb_*

'Test on spillovers

freeze(tb_1) hk3a.wald c(4)=0, c(5)=0,c(6)=0, c(7)=0,c(8)=0, c(9)=0,  c(10)=0
test_bs(1,1)=tb_1(6,2)
test_bs(2,1)=tb_1(6,4)

delete tb_*

freeze(tb_1) in4e.wald c(3)=0, c(4)=0,c(6)=0, c(7)=0, c(8)=0, c(9)=0, c(10)=0, c(11)=0, c(12)=0,c(13)=0, c(14)=0, c(15)=0, c(16)=0,  c(17)=0, c(18)=0
test_bs(1,2)=tb_1(6,2)
test_bs(2,2)=tb_1(6,4)

delete tb_*

freeze(tb_1) ko4ea.wald c(3)=0, c(4)=0,c(5)=0, c(6)=0, c(8)=0, c(9)=0, c(10)=0, c(11)=0, c(12)=0,c(13)=0, c(14)=0, c(15)=0, c(16)=0,  c(17)=0, c(18)=0
test_bs(1,3)=tb_1(6,2)
test_bs(2,3)=tb_1(6,4)

delete tb_*

freeze(tb_1) ma5.wald  c(3)=0, c(4)=0,c(5)=0, c(7)=0, c(8)=0,c(9)=0, c(10)=0,  c(11)=0
test_bs(1,4)=tb_1(6,2)
test_bs(2,4)=tb_1(6,4)

delete tb_*

freeze(tb_1) ph5.wald c(3)=0, c(4)=0,c(5)=0, c(6)=0, c(8)=0,c(9)=0, c(10)=0
test_bs(1,5)=tb_1(6,2)
test_bs(2,5)=tb_1(6,4)

delete tb_*

freeze(tb_1) si3a.wald c(3)=0,c(4)=0, c(5)=0,c(6)=0, c(7)=0, c(9)=0,  c(10)=0
test_bs(1,6)=tb_1(6,2)
test_bs(2,6)=tb_1(6,4)

delete tb_*

freeze(tb_1) ta4aa.wald c(3)=0,c(4)=0, c(5)=0,c(6)=0, c(7)=0,c(8)=0, c(9)=0, c(10)=0,c(11)=0, c(12)=0,c(13)=0, c(14)=0, c(17)=0,c(18)=0
test_bs(1,7)=tb_1(6,2)
test_bs(2,7)=tb_1(6,4)

delete tb_*


freeze(tb_1) th3a.wald c(3)=0,c(4)=0, c(5)=0,c(6)=0, c(7)=0,c(8)=0, c(9)=0
test_bs(1,8)=tb_1(6,2)
test_bs(2,8)=tb_1(6,4)

delete tb_*

